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Abstract. - In this paper, we show that a system of localized particles, satisfying the Fermi 
statistics and subject to finite-range interactions, can be exactly solved in any dimension. In 
fact, in this case it is always possible to find a finite closed set of eigenoperators of the Hamil- 
tonian. Then, the hierarchy of the equations of motion for the Green's functions eventually 
closes and exact expressions for them are obtained in terms of a finite number of parameters. 
For example, the method is applied to the two-state model (equivalent to the spin- 1/2 Ising 
model) and to the three-state model (equivalent to the extended Hubbard model in the ionic 
limit or to the spin-1 Ising model). The models are exactly solved for any dimension d of the 
lattice. The parameters are self-consistently determined in the case of d — 1. 



How do we study systems characterized by the presence of a huge number of interacting 
particles? The standard procedure is a perturbative one. At first, the interaction is completely 
ignored or approximately taken into account as an external field (mean field). Then, a pertur- 
bative expansion is considered. By construction, this scheme works well only when the interac- 
tion is sufficiently weak, but miserably fails when the correlations among the particles are quite 
strong. This is the case for all new highly interacting materials [Fulde(1995)]. An alternative 
procedure has been suggested in Refs. [Mancini and Avella(2003),Mancini(2003)], where the 
main ingredient is the use of composite operators. In the presence of interactions the prop- 
erties of the original (bare) particles are modified and new particles (quasi-particles) should 
be considered. The difficulty with this approach [named composite operator method (COM)] 
is that the number of composite operators rapidly increases with the number of degrees of 
fredom and one must again resort to approximate treatments [Mancini and Avella(2004)]. Re- 
cently, I have found that there is a large class of systems for which it is always possible to find 
a finite closed set of eigenoperators of the Hamiltonian and exact solution can be obtained. 
The purpose of this letter is to define this class, and to illustrate through some examples the 
technique of calculations. Details of calculations and further results will be given elsewhere. 

We consider a system of q species of particles, satisfying Fermi statistics, subject to finite- 
range interactions, localized on the sites of a Bravais lattice. We suppose that the mass of the 
particles is very large and/or the interaction is so strong that the kinetic energy is negligible 
and the particles are frozen on the lattice sites. For a two-body interaction the Hamiltonian 
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takes the form 

H = J2 V * n »® + V M)n a ii)n h (j) (1) 

ia ij ab 

where i stands for the lattice vector and i = (i, t). c a (i) and c\{i) are annihilation and 
creation operators of particles of species a, in the Heisenberg picture, satisfying canonical 
anti-commutation relations. 



{C a (i,*),4(j,t)> = SatSij 

{c a (i,t),c b (j,t)} = {ct(i,t),cl(j,t)} = 



(2) 



V a represents an external field acting on the particles a. n a {i) — c^(i)c a (i) is the charge 
density operator of the fields a. V a b(i,j) describes the two-body interactions, ft is immediate 
to see that the n a (i) satisfies the equation of motion 

i^j) = [ na (i^H} = 0. Then, in order 
to use the equations of motion and Green's function (GF) formalism we must start from the 
Heisenberg equation for the a species 

= VaCail) + ? ^(i, j)n 6 (i)c a (i) (3) 

jo 

where, without loss of generality, we put V aa (i, i) = 0. We see that the dynamics has gen- 
erated another field operator (p a (i) = J2jb Vab(h J) n b(J)c a (i) of higher complexity. By taking 
time derivatives of increasing order, more and more complex operators are generated. These 
operators are composite operators [Mancini and Avella(2003), Mancini and Avella(2004)], as 
they are all expressed in terms of the original fields c a {i) and c\(i). Now, I assert that: be- 
cause of the algebra satisfied by the charge density operator n a (i), the hierarchy 
of composite operators closes. The proof of the statement is the following. For the sake 
of simplicity, let us take Va (i,j) = 2dS a bVaij, where d is the dimensionality of the system, 
ay is the projector on the nearest-neighbor sites and V is the strength of the interaction. 1 
am considering systems with first-nearest neighbor interactions, but the following proof can 
be extended to systems with second, third, ■ • -nearest neighbor interactions. By increasing 
the range of interactions the number of eigenoperators will increase. The time derivative of 
order k is 

c.(t) = V*c a {i) + V mVt m (2dV) m [n a (i)] m c*(i) 



dt 



where n a (i) = ^\ fe aijn (j). Because of ©, the operator [n a (i)] m obeys the following recur- 
rence rule 

2qd 2qd 

K«P = E4 m) KMF E4 ro) = 1 ( 4 ) 

p—1 p—1 

The coefficients A p m ^ are rational numbers that can be easily determined by the algebra 
and the structure of the lattice. Then, for k > 2qd no additional composite operators are 
generated and the equations of motion close. By putting n — 2qd + 1, we construct a n- 
multiplct composite operator 

/ iMO \ 



(5) 



V Mi) J 
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which satisfies the Heisenberg equation 

. dipji) 



l- 



m),H]=ei>{i) (6) 



where the n x n matrix e will be denominated as the energy matrix. Note that we are using a 
vectorial notation: the field ip m (i) is itself a multiplet of rank q. Once the composite operator 
ip(i) and the energy matrix have been determined, an exact solution of the Hamiltonian can 
be obtained. Let us define the retarded Green's function 

G(i,j) = e(t i -t j )({i>(i),^(j)}) (?) 

where (• • •} denotes the quantum-statistical average over the grand canonical ensemble. By 
means of the Heisenberg equation (j^J we obtain in momentum space the equation 

[w-e]G(k,w) = J(k) (8) 

where i~(k) is the Fourier transform of the normalization matrix: = ({ip(i, t), ip'(j, t)}). 

The solution of Eq. © is 

^— 1 ' to - E m + id 

m—l 

where E m are the eigenvalues of the energy matrix e. The spectral density matrices cr( m )(k) 
are calculated by means of the formula [Mancini and Avella(2003)] 



(k)=O om ^O-^/ cb (k) (10) 



where f2 is the nx n matrix whose columns are the eigenvectors of the matrix e. The spectral 
density matrices ^ m \k) satisfy the sum rule ELi £ > (m) ( k ) = M^(k), where Af^(k) 
are the spectral moments defined as 

M^(k) =F.T.({(id/dtr^t),^(i,t)}) (11) 

F.T. stays for the Fourier transform. The correlation function (CF) C(i,j) — \ip{i)ip' (J)/ can 
be immediately calculated from © and one obtains 

n 

C(k, u) = 7T « [w - i? m ]r m a(" 1 ) (k) (12) 

m— 1 

with T TO = 1 + tanh(/3i5 m /2) and /3 = 1/fcgT is the inverse temperature. 

Equations and (|12|l are an exact solution of the model Hamiltonian (TJ- However, the 
knowledge of the GF is not fully achieved yet. The algebra of the field is not canonical: 
as a consequence, the normalization matrix /(k) in the equation of motion (JHJ contains some 
unknown static correlation functions. These correlators are expectation values of operators 
not belonging to the chosen basis and should be self-consistently calculated. 

According to the scheme of calculations proposed by COM [Mancini and Avella(2003), 
Mancini(2003), Mancini and Avella(2004)], one way of calculating the unknown correlators is 
by specifying the representation where the GF is realized. The procedure is the following. 
From the algebra it is possible to derive several relations among the operators. We will call 
algebra constraints (AC) all possible relations among the operators dictated by the algebra. 
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This set of relations valid at microscopic level must be satisfied also at macroscopic level 
(i.e., when expectations values are considered). Use of these considerations leads to some 
self-consistent equations which will be used to fix the unknown correlator appearing in the 
normalization matrix. An immediate set of rules is given by the equation 



1 1 C +oc 



(13) 



where the l.h.s. is fixed by the AC and the boundary conditions compatible with the phase 
under investigation, while in the r.h.s. the correlation function (r/j(i)ip' (i)) is computed by 
means of (|12|) . 

Another set of AC can be derived by observing that there exist some operators, O, which 
project out of the Hamiltonian a reduced part OH = OHq. When Hq and Hj = H — Hn 
commute, the quantum statistical averages of O over the complete Hamiltonian must coincide 
with the average over the reduced Hamiltonian Hq 



Tr{Oe~ 0H } = Tr{0 



} 



(14) 



Another relation comes from the requirement of time translational invariance which leads 
to the condition that the spectral moments, defined by Eq. ifTTjl. must satisfy the following 
relation 

MW(k)=M m W(k)' (15) 

It can be shown that if ljl5|l is violated, states with negative norm appear in the Hilbert space. 
Of course, the above rules are not meant to be exhaustive and, in principle, more conditions 
might be needed. 

We now apply the above procedure to two specific models. It goes without saying that 
the proposed models belong to the class of the Potts models [Potts(1952)]. We consider a 
rf-dimensional cubic Bravais lattice. 

The two-state model 

As first example, we consider only one species of particles. By taking V a = — /i, V a b{i, j) = 
2dV5 a bO-i j the Hamiltonian JQ) takes the form 



H = -n^2v(i) + VdJ2v(i)v a (i 



(16) 



where hereafter for a generic operator <!>(i) we use the notation <E> a (i,i) = J^j aij<i>(j,i). /i is 
the chemical potential and we put v(i) — c^(i)c(i); we are ignoring the subindex a. By means 
of the transformation v(i) = -|[1 + S(i)) the Hamiltonian Ijltjfl can be cast in the form 



H — En 



h"£s(i) 



dJ^S'(i)5 Q (i) 



(17) 



where E = {-\^ + jdV)N, h = |(/i - dV), J = -jdV. Hamiltonian (JT7|l is just the d- 
dimensional spin-1/2 Ising model [Ising(1925)] with nearest neighbor interactions in presence 
of an uniform external magnetic field. 

To solve the Hamiltonian (|16|l let us consider the composite operator 



Mi) 
Mi) 



\ 



( 



\ V>2d+i(i) / 



c(i) 
v a (i)c(i) 



\ 



(18) 
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This field is an eigenoperator of the Hamiltonian (|16|l 

i§- t i>(i) = m),H]=ei>d) (19) 

The energy matrix e has the rank (2d + 1) x (2d + 1) and can be calculated by means of © 
and the recurrence rule @. Proof of the relation @ and the explicit expressions of A$ for 
this model will be given elsewhere. The eigenvalues of the energy matrix are given by 

E m = -fi+(m-l)V {m = l,2, (2d+l)} (20) 

The retarded GF G(i,j) and the CF C(i,j) can be exactly calculated by applying the 
scheme of calculations illustrated above. By using the anti-commutation rules (J5J) and the 
symmetry relations (|15f) . we obtain /(i,j) = 6i_jIo, where the matrix Iq depends on the 
parameters 

K { P ) = (p =!,■■■ 2d). Then, the matrices a (m ^ are calculated by means 

of (|10fl . We have an exact solution of the Ising model for dimensions d = 1,2,3, where all 
the properties of the model are expressed in terms of the self-consistent parameters k^O. To 
have quantitative results we must calculate in terms of the external parameters fi, T and 
V . We now recall the AC (|13|) . Because of the relation Q there are only 2d+ 1 independent 
elements of the matrix C = (ijj(i)ip^ (i)) and we obtain the self-consistent equations 

, 2d+l 

k (p) _ A (p) = _ T m *™ (p = 1, ■ ■ -2d + 1) (21) 

m— 1 

where T m = l+tanh(/3-E m /2) and are the spectral density matrices. New correlation func- 
tions A (p) = (v(i)[v a (i)] p ) appear and the set of self-consistent equations (|21|l is not sufficient 
to determine all unknown parameters. One needs more conditions. In the case d = 1, these 
extra conditions can be obtained by using property 114J) . By means of the algebraic relation 
c^(?)u(i) = 0, we have ^{ije -1311 = c< (i)e~^ H ° , where Ho = H — 2Vv(i)v a (i). By requiring 

that (fUjl be satisfied, we can derive the relation ^ = where G\k = (c(i)c t (i)[w Q (i)] fe_1 ) 

and c[°f} — (c(i)c^(i)[w Q (i)] fc_1 ) . Here (• • -) denotes the thermal average with respect to 
Hq. By exploiting these relations we can arrive to the following equation 

C,3 = i Cl2 (l + gf) (22) 

This equation closes the system of self-consistent equations and we have a complete solution 
of the model. Higher-order correlation functions can be studied by similar techniques. All 
details of the calculations will be given in a separate work, where we will show that all the well- 
known results [Baxter(1982), Goldenfel(1992), Lavis and Bell(1999)] of the one-dimensional 
Ising model are recovered. 
Three-state model 

We consider now a second example. Let us consider two species of particles, say a and 6, 
and take V a — Vb — — /x, V a b(h j) = U8^8 a b + 2dVajj. Then, the Hamiltonian Q becomes 

H = J2[-fin(i) + UD(i) + Vn(i)n a {i)] (23) 

i 

where n(i) — n a (i) + ni,(i) and D(i) — n a {i)nb{i) = ^n(i)[n(i) — 1] are the total particle 
density and double occupancy operators, respectively. This Hamiltonian is just the extended 
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Hubbard model in the ionic limit, where U and V are the on-site and inter-site Coulomb 
interaction, respectively. The two species of particles, a and 6, are in this case electrons with 
spin up and down, respectively. By means of the transformation n(i) = [1 + S'(i)], 1|23[1 can 
be cast in the form 



H = -d.jJ2 S(i)S a (i) + S 2 (i) - hJ^SC*) + E " 



(24) 



where E Q = (-/i + dV)N , h = £t - 2dV -W,J= -dV, A 



\U ■ Hamiltonian 



is just 



the Ising spin-1 model [Blume et al.(1971)] with nearest-neighbor interactions in presence of 
a crystal field A and an external magnetic field h. 

To solve the Hamiltonian (|23() let us consider the composite operators 



( ^ x) W \ 



( x(i) \ 
x(i)[n a (i)} 



\ x{i)\n a {i)} id J 



(25) 



where x = £,?y, being f(i) = [f — n(i))c(i) and r](i) — n{i)c(i) the Hubbard operators. Note 
that c(i), £(i),r](i) are doublet operators with components a and b. By means of (j2J and Q), 
these fields are eigenoperators of the Hamiltonian i|23|) 



= (*)>#] = e ( «V (?) («) 



(26) 



where e^) and are the energy matrices, of rank (4rf+ 1) x (4rf+ f ), which can be calculated 
by means of the equations of motion (|26|) and the recurrence relation The eigenvalues 



-Em'' and £' m , ' ) of the energy matrices are 



(?) 

m 



-\i + (m - l)V 

-fi + U + (m - 1)V 



{m = l,2,---(4d+l)} 



The retarded GF G^(^j') = <-R[^ (a!) (*)^ (w)t 0')]) and the CF ^(i.j) = (V> (a:) W^ (y)t 0')}> 
with x, y = £, 77, can be exactly calculated by applying the scheme of calculations illustrated 
above. By using the anti-commutation relations © and the symmetry relations l)15[l. straight- 
forward calculations show that /^'')(i,j) = /( r '5) (i, j) = , while 



where 
expressions 



(fe = l,..-,4d+l) 

([n Q (i)]P), = §(n(i)[n a (i)p). The CF C^ x \i,j), {x = £,rj), have the 



4rf+l 



(27) 



where Tm = 1 + tanh(/3E™ /2). The spectral density matrices er^™ 1 ) are expressed in terms 
of the elements . We have an exact solution of the model for dimensions d — 1,2,3, 



Ferdinando Mancini: Fermionic systems with charge CORRELATIONS 



7 



where all the properties are expressed in terms of the self-consistent parameters K^and \( p \ 
In order to determine the correlators k^ p ' and \( p \ we use the AC constraints By means 
of the algebraic relations £ c (*)£c W + (*)%(*) = 1 — n c(i), with c = a,b, we obtain for the 
paramagnetic phase the following self-consistent equations 

K (fc-i)_ A ( fc -i) =c[f+c[™ ) (fc = l,.-.4d + l) (28) 

The number of equations is not sufficient to fix all the parameters and more equations are 
needed. In the case of d = 1 these extra conditions can be obtained by using the property 
(1141) . By means of the algebraic relations 

£t(i)„(i) = D p (i)=D(i) 

g(i)D(i)=0 D(i)n p (i)=2 p D(i) { > 



we have 



D(i) e -^ H = D(i) J2 h p [n a (i)]P e - ^ 



where Hq = H — 2Vn(i)n a (i), and the h' p s are known functions of the potential V. By 
requiring that ((Till be satisfied we can derive a set of three self-consistent equations, which 
added to (|28f) allow us to determine all 8 parameters K^and (p = 1, • • -4). Calculations 
and results will be presented elsewhere. 

Summarizing, we have shown that a system of localized particles, satisfying Fermi statis- 
tics, subject to finite-range interactions, can be described in terms of a closed set of eigenoper- 
ators. For the case of nearest neighbor interactions, the number of these composite operators 
is equal to (2qd +1), where q is the number of species of particles and d = 1, 2, 3 is the 
dimensionality of the system. The Green's functions and the correlation functions can be 
exactly calculated and are expressed in terms of a set of self-consistent parameters. For ex- 
ample, we have considered two models: the two-state model (equivalent to the spin-s Ising 
model) and the three-state model (equivalent to the extended ionic Hubbard model and to 
the spin-1 Ising model). For these models the parameters are calculated in the case of d = 1. 
For higher dimensions more self-consistent equations arc needed. This problem is now under 
investigation. 



The author is grateful to Dr. A. Avella for stimulating discussions, a very friendly collab- 
oration and his careful reading of the manuscript. 
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